Observation of Faraday Waves in a Bose-Einstein Condensate 
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Faraday waves in a cigar-shaped Bose-Einstein condensate are created. It is shown that periodi- 
cally modulating the transverse confinement, and thus the nonlinear interactions in the BEC, excites 
small amplitude longitudinal oscillations through a parametric resonance. It is also demonstrated 
that even without the presence of a continuous drive, an initial transverse breathing mode excita- 
tion of the condensate leads to spontaneous pattern formation in the longitudinal direction. Finally, 
the effects of strongly driving the transverse breathing mode with large amplitude are investigated. 
In this case, impact-oscillator behavior and intriguing nonlinear dynamics, including the gradual 
emergence of multiple longitudinal modes, are observed. 

PACS numbers: 03.75.Kk,67.90.+z,05.45.-a,47.54.-r 



In 1831, Faraday studied the behavior of liquids that 
are contained in a vessel subjected to oscillatory verti- 
cal motion [l[. He found that fluids including alcohol, 
white of egg, ink and milk produce regular striations on 
their surface. These striations oscillate at half the driv- 
ing frequency and are termed Faraday waves. They are 
considered to be an important discovery. Since then, the 
more general topic of pattern formation in driven systems 
has been met with great interest, and patterns have been 
observed in hydrodynamic systems, nonlinear optics, os- 
cillatory chemical reactions, and biological media [2]. 

In this paper we study pattern formation by mod- 
ulating the nonlinearity in a Bose-Einstein condensate 
(BEC). Nonlinear dynamics arise from the interatomic 
interactions in this ultracold gas. In the past the ob- 
servation of interesting phenomena has motivated re- 
searchers to propose and implement various techniques 
to manipulate the nonlinearity. Such control has been 
accomplished for example by exploiting Feshbach reso- 
nances In our experiment we investigate an alterna- 
tive technique, namely periodically modulating the non- 
linearity by changing the radial confinement of an elon- 
gated, cigar-shaped BEC held in a magnetic trap. The 
radial modulation leads to a periodic change of the den- 
sity of the cloud in time, which is equivalent to a change 
of the nonlinear interactions and the speed of sound. This 
can, in turn, lead to the parametric excitation of longitu- 
dinal sound-like waves in the direction of weak confine- 
ment. This process is analogous to Faraday's experiment 
where the vertical motion of the vessel produced patterns 
that were laterally spread out. 

It has been shown theoretically that for a BEC, a Fara- 
day type modulation scheme in the case of small driv- 
ing frequencies leads one to the same type of analysis as 
would the direct modulation of the interatomic interac- 
tion, e.g., by a Feshbach resonance 

US]. In both cases, 
the dynamics are governed by a Mathieu equation that is 



typical for parametrically driven systems. Floquet anal- 
ysis reveals that a series of resonances exist, consisting 
of a main resonance at half the driving frequency, and 
higher resonance tongues at integer multiples of half the 
driving frequency 

In our experiment we exploit this transverse modula- 
tion scheme for three different applications. First, we 
apply a relatively weak continuous modulation, demon- 
strate the emergence of longitudinal Faraday waves, and 
study their behavior as a function of the excitation fre- 
quency. Second, we investigate longitudinal patterns that 
emerge as a consequence of an initial transverse breath- 
ing mode excitation without the presence of a continuous 
drive. This has important consequences in the context of 
damped BEC oscillations and has been studied theoret- 
ically in p. Since the first experiments with BECs, the 
study of collective excitations has been a central theme 
[7|, [8[. The transverse breathing mode, which we exploit 
in our experiments, plays a prominent role: Chevy et al. 
[9[ showed that this mode exhibits unusual properties, 
namely an extremely high quality factor and a frequency 
nearly independent of temperature. Finally, in a third 
set of measurements we study the situation of a rela- 
tively strong modulation, resonantly driving the trans- 
verse breathing mode. We show that the condensate re- 
sponds as an impact oscillator, which leads to intriguing 
multimode dynamics. 

The experiments were carried out in a newly con- 
structed BEC machine that produces cigar-shaped con- 
densates of 87 Rb atoms in the \F = l^mp — — 1) 
state. The typical atom number in the BEC is 5 • 10 5 , 
and the atoms are evaporatively cooled until no thermal 
cloud surrounding the condensate is visible any more. 
The atoms are held in a cylindrically symmetric Ioffe- 
Pritchard type magnetic trap with the harmonic trap- 
ping frequencies of {(j xy /(27r),(j z /(27r)} = {160.5, 7}Hz. 
The weakly confined z-direction is oriented horizontally. 
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For the experiments described below, the following col- 
lective mode frequencies are of particular importance: 
first, there exists a high-frequency transverse breathing 
mode. For our trap geometry, this mode has a frequency 
of co±/(2tt) = 321 Hz [10|, [ill], very close to the limit of 
vanishing axial confinement 2 • uo xy / (2tt). The second set 
of modes in which we are interested here consists of axial 
modes which, for large quantum numbers, correspond to 
sound waves in the z-direction. The frequencies of this 
discrete set of modes can be approximately calculated as 
given in TH, 13j |. 

In order to investigate the parametric driving pro- 
cess mentioned above, we first performed a set of "spec- 
troscopy" experiments in which we continuously mod- 
ulated the transverse trapping confinement at a fixed 
modulation frequency and observed the subsequent emer- 
gence of longitudinal Faraday waves. For each excitation 
frequency the modulation amplitude was adjusted such 
that the longitudinal patterns emerged typically at some 
point after 10 to 30 oscillations. On the breathing mode 
resonance, a trap modulation of 3.6% was used, while at 
many other frequencies trap modulations of up to 42.5% 
were chosen to obtain clearly visible patterns [14(. This 
range of modulation depths is similar to the range used 
in numerical simulations in [5[. Representative exam- 
ples for the resulting patterns are shown in Fig. [TJ All 
experimental images in this manuscript were taken by 
destructive in-trap imaging. 

The average spacing of adjacent maxima in the result- 
ing pattern is plotted against the driving frequency in 
Fig. [21 The data lie on a clear curve, with the exception 
of the points near a driving frequency of 160.5 Hz, cor- 
responding to the transverse dipole mode resonance (i.e. 
transverse slosh motion) [To| . However, inspection of our 
experimental images reveals that, at this frequency, we 
also excite the transverse breathing mode at 321 Hz. Ex- 
citation of the breathing mode is a very effective way of 
creating longitudinal patterns. Therefore the patterns 
obtained at 160.5 Hz are actually the same as those pro- 
duced at 321 Hz. In order to rationalize the data, we first 
note that parametric excitation with a certain driving 
frequency excites oscillations predominantly at half the 
driving frequency, the main resonance also observed in 
Faraday's experiments. The dispersion relation of longi- 
tudinal collective modes that become sound-like for high 
quantum numbers is given in [12|, |l3[ and is used to calcu- 
late the expected spacing between density maxima. The 
resulting spacings are plotted as the step-like curve in 
Fig. [2] and are in excellent agreement with our experi- 
mental data, corroborating the assumption and theory 
of a parametric driving process. 

In a second set of experiments we show that the lon- 
gitudinal modes are driven by the transverse breathing 
motion even without the presence of a continuous exter- 
nal drive. In particular, this disproves the influence of 
any tiny residual axial trap modulation that has existed 
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FIG. 1: In-trap absorption images of Faraday waves in a 
BEC. Frequency labels for each image represent the driving 
frequency at which the transverse trap confinement is modu- 
lated. 
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FIG. 2: Average spacing of adjacent maxima of the longitudi- 
nal patterns plotted versus the transverse driving frequency. 
Points are experimental data, while the line shows the theo- 
retical values calculated for the longitudinal modes closest to 
half the driving frequency. 



during the presence of our continuous drive due to exper- 
imental imperfections [l4(. For this, we excited the trans- 
verse breathing mode at 321 Hz by driving the transverse 
trap confinement for a few cycles, and then let the con- 
densate evolve without the presence of the drive. For 
a gentle excitation, we do not see longitudinal patterns 
immediately after the end of the drive, but can observe 
them emerge at later times. A weaker excitation delays 
the emergence of longitudinal modes out to later times, 
while in the case of strong excitations the patterns can 
emerge within the first three cycles. In order to follow 
the evolution, we quantify the presence of longitudinal 
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FIG. 3: Onset of Faraday wave formation after driving the 
transverse breathing mode for only two cycles. Pattern visi- 
bility is quantified as described in [ToT ]. 

patterns as described in [15|. Fig. [3] shows the evolution 
after the end of a moderate excitation. For this data, we 
excited the condensate for only two cycles, varying the 
transverse trap frequency by 9% with a modulation fre- 
quency of 321 Hz during the modulation. Immediately 
after the excitation, the obtained images showed no longi- 
tudinal waves, and our pattern visibility measure, plotted 
in the figure, is initially picking up high frequency noise 
along the longitudinal axis. A main contribution to this 
noise is the imaging noise of our detection system. Weak 
pattern formation is observed starting at five periods af- 
ter the end of the modulation, and strong longitudinal 
patterns then appear after about nine periods. A similar 
behavior is known for example from parametric ampli- 
fiers in optics: if no input signal amplitude is present, 
a signal emerging from noise (or zero point energy) can 
form if the amplification is large enough. This behavior 
is also found by our numerical simulations based on the 
Gross-Pitaevskii equation. In simulations with similar 
parameters as in the experiment, no pattern formation is 
observed for several breathing mode periods. From the 
onset of pattern formation, it takes just three periods for 
patterns to grow to their full strength. We find that the 
onset time of pattern formation in the simulations is ear- 
lier when more noise is added to the initial relaxed wave 
function. 

This experiment is closely related to the theoretical sit- 
uation described in where a single and sudden jump 
in the transverse trap frequency was used to excite the 
breathing mode, instead of a sinusoidal trap frequency 
modulation. We have used a trap jumping excitation 
in the case of very elongated BECs to demonstrate a 
second observation about the onset of longitudinal pat- 
terns, namely the fact that these patterns can start to 
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FIG. 4: Pattern domains observed in a very elongated BEC. 
The graph shows an integrated cross section of the BEC dis- 
played in the inset. 

emerge in spatially localized domains, rather than uni- 
formly across the whole cloud. To show this effect, we 
produced condensates in a very elongated cigar-shaped 
trap with trapping frequencies of {u xy /(27r), u z /(2tt)} = 
{286.1, 2. 8}Hz. We temporarily jumped to a different 
trap of {u xy /(27r),u z /(27r)} = {88.4,5.1}Hz for the du- 
ration of 1.3 ms, then jumped back to the first trap and 
let the cloud evolve. After about 10 ms, we observed 
BECs in which a perfectly periodic density modulation 
was stretching over almost the entire cloud; but it is also 
not uncommon to observe patterns in several separate 
domains, as shown in Fig. [U Considering that the speed 
of sound in the elongated cloud is only about 2 mm/s 
[l3L [l6| ], it is plausible that over the time scale of this 
experiment different regions of the BEC can get inde- 
pendently excited and evolve into longitudinal patterns 
independently from each other. 

In the experiments described so far, our parametric ex- 
citation has led to longitudinal modes oscillating at half 
the frequency of the parametric driving frequency. In 
the case of strong parametric amplification, it is theoret- 
ically expected that modes at other frequencies (higher 
resonance tongues), in particular modes at the driv- 
ing frequency, can be excited, too [4]. This motivated 
our third set of experiments in which we started again 
with a condensate in a trap with trap frequencies of 
{cjs2//(27r), 0^/(271-)} = {160.5, 7}Hz. We then continu- 
ously modulated our transverse trap frequency by about 
19% with a modulation frequency of 321 Hz. The result- 
ing breathing motion is seen in Fig. [5] where we plot the 
Thomas- Fermi radius of the cloud in the transverse direc- 
tion versus time. The graph clearly shows that the cloud, 
upon strong excitation, gradually starts behaving as an 
impact oscillator, i.e., as an oscillator bouncing off a stiff 
wall during each period. A classical impact oscillator, re- 
alized for example by a ball bouncing off a stiff surface, is 
a paradigm for nonlinear dynamics, nonlinear resonances 
and chaotic behavior. In the present case, the role of the 
stiff wall was played by the strong mean-field repulsion 
during the slim phase of the oscillations when the BEC is 
strongly compressed in the transverse direction. Theoret- 
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FIG. 5: Impact oscillator behavior. Plot of BEC radius in 
the tight direction vs. time when continuously driving the 
transverse breathing mode. Points show experimental data, 
while the curve corresponds to numerical simulation. Inset: 
(a) Image of BEC after 5.2 driving periods, i.e. 16.2 ms. (b) 
Fourier transform of the image in (a). 



ically, the instability of the breathing mode upon strong 
driving and the impact oscillator behavior have been an- 
alyzed in [17]. The behavior of the cloud radius in the 
transverse direction was also reproduced in our numer- 
ical simulation of the azimuthally symmetric 3D Gross- 
Pitaevskii equation, displayed as the solid line in Fig. [5j 
However, the numerics, when starting with a thoroughly 
relaxed wave function in the initial trap, show a sign of 
longitudinal pattern formation only after 18 ms, while in 
the experiment, longitudinal patterns clearly formed al- 
ready during the third period (9 ms). This, again, hints 
at the importance of initial noise in the condensate that 
seeds the parametric amplification. In the experiment, 
the patterns start out similar to those displayed in Fig. [1] 
for the case of a weak drive at 321 Hz. But, upon the 
action of the strong drive, they quickly evolve into more 
complicated patterns, involving the excitation of several 
other modes. The inset of Fig. [5] shows an image taken 
after 5.2 driving periods (a), together with the Fourier 
transform (b). The Fourier spectrum reveals that several 
modes corresponding to the first resonance tongue of lon- 
gitudinal modes with nearly half the driving frequency 
are excited. In addition, modes at twice the distance 
from the central Fourier peak are visible. Those modes 
belong to the second resonance tongue of the main reso- 
nance. 

In conclusion, we have experimentally observed the ef- 
fects of parametric resonances in a BEC. The observed 



resonances lead to Faraday waves along the long BEC 
axis. These results advance the understanding of col- 
lective mode behavior in a condensate, which is one of 
the key tools to study BEC dynamics. In addition, we 
have shown that strongly driving the transverse breath- 
ing mode leads to an instability whereupon the mode 
amplitude increases exponentially, accompanied by the 
strong excitation of multiple sound-like modes. 
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